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Jesse Pino. 1 '!*^ Hui Li, 1 and Swadesh Mahajan 2 

1 ' Los Alamos National Laboratory, Los Alamos, NM 87545 

2 ' Institute for Fusion Studies, University of Texas at Austin, Austin, TX 78712 
(Dated: 4 November 2010) 

The evolution equations for a plasma comprising multiple species of charged fluids with relativistic bulk and 
thermal motion are derived. It is shown that a minimal fluid coupling model allows a natural casting of the 
evolution equations in terms of generalized vorticity which treats the fluid motion and electromagnetic fields 
equally. Equilibria can be found using a variational principle based on minimizing the total enstrophy subject 
to energy and helicity constraints. A subset of these equilibria correspond to minimum energy. The equations 
for these states are presented with example solutions showing the structure of the relaxed states. 
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I. INTRODUCTION 

In a variety of astrophysical settings such as jets pow- 
ered by Gamma Ray Bursts, supermassive black holes, 
or pulsars, energies can be such that the bulk motion of 
a fluid approaches the speed of light, and/or the ther- 
mal motion is comparable to the rest mass of the par- 
ticles that make up that fluid. These cases necessi- 
tate a relativistic treatment of plasma motion. Previ- 
ously, Mahajan- showed that the dynamics of relativistic 
charged fluids could be treated using a minimal coupling 
formalism in direct analogy to the canonical momentum 
prescription of single particle dynamics. This formal- 
ism has been successfully used to describe both linear 
and nonlinear waves in relativistic pair plasmas of pul- 
sar magnetospheres 2 ^. Presently, we demonstrate that 
there is a natural minimization principle associated with 
the minimal coupling formalism. When total energy and 
the helicity of multiple charged species are conserved, 
the application of this minimization principle to config- 
urations leads to relaxed equilibrium states. We believe 
that these states give a better description of the physical 
system and can be used as more accurate starting points 
for subsequent dissipation and slow evolution. 

The outline of this paper is as follows. In Section 2, 
we present the formalism behind relativistic magnctoflu- 
ids and derive a minimization principle for finding equi- 
librium states. In Section 3, we present some examples 
of one dimensional solutions to the equilibrium equations 
derived from the minimization principle. Finally, we con- 
clude by summarizing our results and discussing the ap- 
plicability of these relaxed states to astrophysical phe- 
nomena, namely the 'striped wind' of a pulsar nebula 4 . 



II. RELATIVISTIC CHARGED FLUIDS 

We consider an isolated system of multiple fluids, and 
take the mass-energy to be small enough that we can ne- 
glect changes to the geometry of space-time. We thus 
take the Minkowski metric g^ v — diag(—l, 1, 1, 1). The 
velocity 4- vector is = (7c, 711), where u is the local 3- 



velocity of the each fluid species and 7 = (1 - 

„2 



7c 2 ) 



-1/2 



Then U 2 = U^U^ = —c z . We assume local Maxwellian 
closure, so all of the fluid quantities are implicitly func- 
tions of position (e.g. n, T, p), with n the proper density 
of the fluid, and the pressure p = nkT. 

Each species (labeled s) satisfies mass conservation, 

d v V v s = d v {n s U v s ) = 0, (1) 

and has the stress-energy tensor 

Tr = p s g^ + n s m s G{z s )U^U v , (2) 

where nmc 2 G{z) = p + ne is the enthalpy (pressure plus 
internal energy) density of the fluid£, and we assume 
some fluid closure that gives the enthalpy as a function 
temperature (z — mc 2 /kT). The results presented in 
this section are independent of the specific form of this 
enthalpy. Each fluid has a charge q s , coupling it to the 
electromagnetic field through the Lorentz force equation, 



(3) 



where the E.M. field tensor F» v = d^A v ~ d v A» , with 
A^ = (<f>, A) the 4-potential 6 . The E.M. field also obeys 
Maxwell's equations (F is the dual of F), 



4?r 



J" , d v F^ 



0, 



and has stress-energy tensor 
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— ( n^ a F aF^ v + - a ^"F a ^F , 



(4) 
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In order to close the system, we define the current 
J v = ^2 s q s T" as the sum of the fluid currents, which 
is required so that the total stress-energy T^ t = T^ M + 
£ s is divergence free {d v T% t =0). 

Expanding the force equation ([3]), we obtain (we sup- 
press species labels except when needed for clarity): 



1 



U v d u (mGU^) + -d" P = qF^U„ 



(5) 



Using U u d v = d/dr = jd/dt = j(d/dt + u • V), we can 
break this up into components: 



— (mc^G) — p = qu • E. 

dt n-f at 



(6) 



4(mG7u) + — Vp = g(E + ux B/c). (7) 

dt 717 

This is the relativistic generalization of the Lorentz force 
equation, and we see that the effective mass is increased 
both by bulk motion (7) and thermal motion (G)2r— . 
If we define the entropy a through 



Td v a = mc 2 d u G--d"p, 
n 



(8) 



we can rewrite the evolution equation ((3J) for each species 
as 

qU„(F^ + (c/q)S" v ) = qU„M^ = Td v o (9) 

where S» v = d^{mGU u ) - d u (mGU^) is an antisym- 
metric tensor constructed using the temperature trans- 
formed momentum as a potential - 1 - ^" . This is the 'Mini- 
mal Coupling Magnetofluid Unification' model described 
in Mahajan (2003) 1 . The justification is that the com- 
bined potential W = A v + (mc/ q)GU v can be seen as 
the canonical momentum of the fluid, just as one has 
p — s- p + qA for single particle motion. Since M^ v is 
anti-symmetric, contracting equation (J9j> with gives 
da I dt = 0, the standard result of constant entropy along 
field lines. 

By using fluid potential mGU 1 ' = (mcGj, rnGju) = 
(x, P) in place of <fi and A, we can define the fields Q — 
(—dtP/c — V%) and R = V x P as direct analogues to 
the electromagnetic fields E and B. This allows us to 
expand © as 



qu- E 



c \ _ Tda_ 
q J j dt 



(10) 



4>+(mc 2 /q)G'y. Equation (p~0|) states that entropy is gen- 
erated locally only if the fluid flow is not perpendicular to 
the sum of the electric field and the inertial analogue of 
the electric field. Equation (TTTT) describes the evolution of 
the magnetic field and fluid momentum, including both 
gradient and inductive forces. Taking its curl, we ob- 
tain the generalized vorticity evolution equation for each 
species: 



dt 



v x (u s x n s 



X VfT. 



(12) 



The first two terms in equation (| 1 2[) state that the gen- 
eralized vorticity is 'frozen-in' to the flow. The last 
term describes vorticity generation due to the relativis- 
tic extension of the baroclinic term. To make con- 
tact with the familiar MHD equations, we take a mass- 
less barotropic electron fluid with u as the bulk plasma 
flow. We then recover the ideal MHD induction equa- 
tion <9 t B = Vx(ux B), and the entropy generation rate 
Td t o = J E . 



III. MINIMIZATION PRINCIPLE 

We now seek equilibrium configurations of systems 
with multiple species. The presence of dissipation will 
allow ideally conserved quantities to change over time. 
In general, the conserved with the highest order spatial 
derivatives will change the fastest, and those with lower 
order derivatives can then be treated as constraints in a 
variational principle with the most susceptible quantity 
as the target of minimization 12 . A standard example 
in MHD is the minimization of magnetic energy (B 2 ), 
subject to the constraint of constant magnetic helicity 
(A • B ~ B 2 /k): 



8{E-ah) = 0. 



(13) 



This results in the well-known Taylor state V x B = 
aBi£. Within the current framework, species helicity 
and total energy are both ideally conserved variables un- 
der suitable assumptions (see Appendix). However, the 
kinetic terms cause the total helicity to include higher- 
order spatial derivatives than the total energy, making it 
more fragile to dissipation. Thus the variational princi- 
ple (|13l) is mathematically ill-posecU^. Simply minimiz- 
ing helicity subject to constant energy is not possible ei- 
ther, since in general the helicity is not positive-definite 
and often not bounded from below. Instead, we must 
find a coercive functional N which is fragile to dissipa- 
tion but has a quadratic form. We can then minimize N 
while keeping both species helicity and total energy as 
constraints: 



0_ 

dt 



T 

uxJ7-V$H Vcr, (11) 
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S JV - fiE -J2 a s h * I = °- 



(14) 



where we have defined the total species vorticity ft = 
V x [A + (c/q)P] = B + (c/q)R, and potential $ = 



The resulting states can then be compared to the equilib- 
rium conditions to make sure they are physical, and we 
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can further seek a subset of these relaxed states which 
correspond to minimum energy. 

We begin by constructing a canonical stress tensor for 
each species by replacing F in equation Q with the 
species canonical tensor M^ v = F^ v + (c/q)S^ v : 



4ir 



(15) 



Assuming that the system is isolated and that the fields 
vanish at the boundaries, the integral of the time-time 
component of this tensor will be an effective Hamiltonian 
for this system^. Summing over species, we see that this 
is equivalent to the sum of the species total (magnetic 
plus fluid) enstrophies 



7V = V- 



Sit 



|V$| 



la 



d 3 3 



(16) 



The enstrophy is the most fragile to dissipation, since 
it contains the square of the curl of the fluid momentum 
(R terms in ft). This makes it a proper target functional 
for minimization. We will minimize this with constraints 
of the total energy and species helicities, by varying the 
potentials <f> and A and the momentum P s . Note that 
we do not independently vary x, the time component of 
GU^, since at fixed temperature, 



S X = 5(cGj) = S{y / c 2 G 2 + P 2 ) 
The variation of the enstrophy is: 



SP 



SP. 



[V x Sl s ] ■ SA 



+- 



m s c 



Vxfl, 



■5P S }d 3 x.(17) 



The total energy is the sum of electromagnetic and 
kinetic energies 



E = I T 0Q d 3 x = E EM + E, 
E 2 +B 2 



(18) 



Stt 



-d 3 x 



J (mnc 2 j 2 -p) s d 3 



which leads to the variations: 



SE EM = /^{(-V 2 W+(VxB)^A}A, (19) 



SE S — SI n s m s c z G s 7 2 d 3 x 



n,in»S ( ° Gs ^~ P<i ) <f'.r 



= I 2n s m s j s u s ■ SP S d x, 



(20) 



We define the generalized helicity in the Appendix, as 

h = J[B + (c/q)R]-[A+(c/q)P}d 3 x, (21) 

which contains the fluid helicity, magnetic helicity and 
cross helicity terms. The variation of this helicity gives 



6h. 



n s ■ SA + —n s ■ sp. 



(22) 



Plugging these variations into the minimization equa- 
tion (|14p and considering each variation independently 
results in the set of equations: 



^V 2 $ s -/iV 2 = O, (23) 



E 



47T 



-V x fL — a,fi s 



47T 



V x B = 0, (24) 



4tt L 



u 



Vxfl, -V 2 $ 



2(j,n s q s j s — (25) 
c 

—a.tl. = 0. 



To ensure that these solutions are physical, (|2"31 and (f2~4")l 
must be equivalent to Poisson's and Ampere's equations, 
which requires 



Y V 2 $ s = -47TfiJ2n s q s 



Substituting ([53J into (J23) 
1 



JL 

-in 



4tt 



V 2 $, 



2un s q s 



Us 

c 



(26) 



(27) 



Thus each species component of (|26[) must hold indepen- 
dently, viz. 



1 „ 9 / , m s 

-v 2 U + — Xl 



-4:Trn s q s j s . 



(28) 



A* V Is 

The states with minimum enstrophy satisfy the relation 

1 „ „ n,o,7, 

-Vxft.-aA^/i^u, 29 

47T C 

The subset of minimum energy states can be found by 
taking the limit fx, a s oo, keeping a s /fi = — 1/A S con- 
stant. This formally reproduces the energy minimization 
S (E — J2 S h-s/^s) — 0, however, by using the states con- 
structed from (fl4)) . we arrive at the relaxed states via 
a well-posed minimization problem. In order for both 
Ampere's law and (|2"8|) to still be satisfied in the limit of 
large \i, we must have 



0. 



(30) 



For a system with V(</> + ^Xs) = on the boundaries, 
this implies that 



7s G s 



m s c* 



t4> = In 



const, 



(31) 
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which means that the gradient forces in the evolution 
equation vanish. Since G, 7 > 1, mc 2 ^ max can be in- 
terpreted as the total energy available to the fluid species, 
to be distributed among thermal, bulk motion, and elec- 
trostatic energies. 

The relaxed fluid states further satisfy the relation 



V x (G 7 u) s = A S ^4V, 



m s cr 



-^B. (32) 



m s c 



This condition was given previously by Elsasser and 
Pope! 1 '' in the case of non-relativistic temperatures, with- 
out justification for it being a minimum energy state. 
This equilibrium condition introduces a natural length 
scale of the collisionless skin depth. To show this more 
concretely, we take the case of an electron-positron pair 
plasma. Normalizing length to some system size L, ve- 
locity to c, density to no, and magnetic field to Bq = 
V27rm e 7i c 2 , Equations ([27]) and (|32|) become 



eV x P , 



eV x B 





P+ 


B 

~ 2"' 


(33) 


6 A "G_ 


P 


B 

+ r 


(34) 


n+ P 


n_ 

gZ 


p , 


(35) 



where A = X/8wL, and e = X pe /L = ^/ mc 2 / 8irnoq 2 / L 
is the ratio of the pair skin depth to the system size. 
P + = G + 7 + u + is the momentum of the positrons, and 
minus sign subscripts denote electron quantities. Previ- 
ous work by Iqbal et alJ^ examined the limit of isotropic 
equal temperatures and densities, with quasineutrality 
and nonrelativistic bulk flow (7 ~ 1). They showed that 
B could be expressed as the sum of three Beltrami flows 
(V x Fj = fJ-jFj), with different scale lengths ^J 1 . In ap- 
propriate limits, one of the eigenvalues fi can become very 
small, leading to structures on a scale much larger than 
the skin depth, making the solutions relevant to mod- 
elling of astrophysical jets. Indeed, if we assume that the 
spatial variation of the fields is only on the large scale, 
to first order we recover the MHD relaxed "Taylor state" 
V x B = aB, with a = (A^ 1 + Al 1 )/2. However, the 
presence of the small parameter e multiplying the curl 
operator means that the multi-fluid equations represent 
a singular perturbation on MHD dynamics, and that the 
small scale structure cannot be fully ignored. In general, 
the relations between G, 7 and n make the equilibrium 
equations (|33H35I) highly nonlinear, and solutions must 
be found numerically. In the following section, we ex- 
amine the class of one dimensional solutions in various 
simplifying cases. 



IV. ONE DIMENSIONAL SOLUTIONS 

In order to explore some configurations of the re- 
laxed states, we will now restrict to the class of solu- 



tions in one dimension. We take a Cartesian coordi- 
nate system with all quantities varying in x only. We 
further consider those simple cases where the enthalpy, 
temperature and density of the two species are equal 
(G+ = G_ = G, T+ = T_ = T, n+ = n_ = n). This 
in turn implies that A+ = \ 2 _ (We drop the hat notation 
of the previous section). For definiteness, we take the 
relativistic Maxwellian distribution ^ 18 i 19 . The enthalpy 
then takes the form G(z) = Ks{z) / l Ki(z) , where Ki is the 
MacDonald function of order i and z = mc 2 /ksT. We 
also take the entropy to be constant within the region of 
interest so the density takes the form n = uq (K2 (z)/z). 
This assumption can later be relaxed (see the following 
discussion). 



A. Perpendicular current 

We first look for an analogue of the 1-D Harris-Hoh 
shee t 20 ! 21 , with JIB. This implies that the direction of 
B does not change, and we are free to choose B = B(x)z, 
and J = J{x)y — —B'y. We could take in addition a con- 
stant J x o, but then the current would not vanish at the 
boundaries. Also note that the 1-D assumption precludes 
the possibility of a constant B x q, since equations (|33[) and 
(1341) would demand nonzero curl of P in the x-direction. 
The current equation (|35|) then becomes 



dB n . 

e a^-G (p + 



p-), 



(36) 



implying that P x _ = P x+ , and P z _ = P z+ . 
Crossing (|33|) and 134[) with B, we obtain 



eV(P + ■ B) = e(B • V)P+ + e(P+ ■ V)B 



- e -i^-(BxP H 



eP+x(VxB), (37) 



eV(P_ ■ B) = e(B ■ V)P_ + e(P_ ■ V)B 
-e -1 A_5(B x P_) + eV x (V x B). 



(38) 



The first term on the right implies that P x = for both 
species. Then by inspection, the above equations are 
equivalent only if we take A+ = — A_ = A, and P+ = 
P,| + P x = P z z + P y y, and P_ = P,, P ± = P z z - P y y. 
Thus there is a net current perpendicular to the magnetic 
field as well as net momentum parallel to the field. The 
current equation is then written 



e ^--2^P 
6 dx ~ G 



Equations (j3"Tf and (|3"8")l become 



2P Z )P 



11 



(39) 



(40) 



Using (|39|) allows us to write: 
,dP z 



eB- 



dx 



2-P P 
G 



v 



G 



(e- x \B + 2P Z ) P y 
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or, 



dP z 
dx 



(41) 



Dotting fl33j) and (J34J) with B gives 
eV ■ (P ± x B) = 



2 



2 G P 



2 

±5 



which in 1-D reduces to 
Using (J3SJ), we get 

<9x 



2 



I? 
2~ 



2—P 2 

G v 



(42) 



We have a set of three equations (J22]) , dH]) , and (g2]) . 
In order to close these equations, we take the adiabatic 
distribution where 7G = "f m ax- Then we are able to 
write G 2 = 7^j OX — P 2 ■ The temperature and density can 
be solved for using the relativistic Maxwellian distribu- 
tion described above. The nonlinearity inherent in this 
prescription necessitates numerical integration to find so- 
lutions. As an example, in Figure [I] we plot a 1-D equi- 
librium solution for a Maxwellian pair plasma containing 
a magnetic field reversal similar to a Harris sheet. The 







^tot 




T 








J y 






X 



FIG. 1. Example of a pair-plasma equilibrium state in a 1-D 
B z field-reversal configuration, with A+ = — A_ = 1. Per- 
pendicular to the field, the electrons and positrons stream in 
opposite directions, generating a current J y . Along the field, 
both species flow together, creating a parallel flux F z = wyu z . 
The adiabatic condition leads to peaked density (ntot = 
n+ + n_) and temperature profiles, with points of vanishing 
density, beyond which the magnetic field is constant. Distance 
is normalized to the skin depth, and n — no, fcsT = mc 2 at 
x = 0. 

electrons and positrons have equal temperatures and den- 
sities, but opposite eigenvalues A + = — A_ = 1. The adi- 
abatic condition pip leads to a compact area of nonzero 
density and current, outside of which the magnetic field 



is constant. The finite helicity of the system allows for 
bulk flow along the magnetic field. The two-fluid nature 
of the equations leads to field reversal on the scale of the 
skin depth, suggesting application to magnetic reconnec- 
tion studies and shock structure. The accessibility and 
stability of such nonlinear multi-fluid equilibria will be 
explored in forthcoming work. 



1. Oscillatory and Constant Solutions 

Differentiating eq (l4"2"l) . we have 



e^AQ)'p, + e -^-4 



e 2 P" = A 



(5)' 



P Z ~ 



-2 A 2 



- 1 



G 



Py (43) 



If the variation in n/G is small, then the solution is os- 
cillatory when 



A > ey/G/n = 



1 



C /Up,eff 



L y/ 'STrmnofiq 2 



where w P , e // is the effective plasma frequency including 
thermal effec ts. 

If A = e^J G/n, and P v = 0, there is a solution with 
constant B,P z ,n, and T. From equation (|4"2")l , we see 
this requires P z = \J G/nB/2, or the total flow velocity 
(electrons plus positrons) is u z = B/y/nGj 2 . That is, 
the magnetic field is constant with flow along field lines 
at a relativistically modified Alfven speed. 



2. Pressure Balance 

Examining Equations ([4"Tj) and (|4"2"j) . we see that 
we have 



dB 2 n dP 2 
= 4- 



dx 



G dx 



(44) 



If we take the adiabatic distribution 7G = 7 ma :r, we are 
able to write G 2 = j^ iax — P 2 . At x = 0, we take B = Bq, 
T = Tq, and n — uq. The boundary xx is then defined 
as the point where the density vanishes and P = P m ax = 
\Jlmax 1- For x > x\, the magnetic field is constant, 
B = B(xi) = B\. Then equation (|44"]) implies that 



B 2 



B 2 



dB 2 

dx 



dx 



(45) 



ndP 2 , 
i G^x~ dx = 



n BP 2 dG 
G 8G dx 



Using the fact dP 2 /dG — — 2G, and that the constant 
entropy distribution has 

dG _ 1 d{nT) 
dx n dx ' 
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we arrive at (since n\ = 0) 



B\ - B 2 = 



djnT) 
dx 



dx = 8n.oTo. 



(46) 



Thus the difference between the boundary magnetic field 
and the central magnetic field strength depends only on 
the central density and temperature, but is independent 
of A, the width of transition, or the orientation of the 
fields. This is equivalent to saying that the magnetic 
pressure and the thermal pressure balance. 



B. Equal A: no bulk motion 

If A+ = A_ = A, we can have a solution with no total 
momentum, Vtot — P+ — P- = 0, and equations (|33|) 
and (|34p become degenerate. Thus 



eVxP = e-^P-|, 
where P = P + = P . In 1-D. this reduces to 



dP z 
dx 



, n 1 



dx " "G 
The current equation (|35|) becomes 
dB 



dx 



2-P 
G 



dB, 
dx 



— —2—P 

- Q V 



(47) 



(48) 



(49) 



(50) 



(51) 



We note that the 1-D constraint enforces P x — B x = 0. 

In Figure [2j we show two examples of one dimen- 
sional solutions with equal A = 1. Both solutions have 
To = 1, no = 1, and P y (xo) = 1. In the first plot (a), the 
central magnetic field vanishes at the origin. There is a 
finite region over which the plasma density is supported. 
The magnetic field changes in both direction and magni- 
tude over a width of Ax = 9.77. As discussed above, the 
magnitude of B at the edge is set by the central density 
and temperature. In this case, B 2 (x\) — 8.0, and the 
field orientation twists through an angle of A6 = 131.5°. 
In (b), the central field is taken as B y (xo) — 1. Nowhere 
does the magnetic field vanish in this case. The transi- 
tion width is Ax = 8.68, the edge magnetic field strength 
is B 2 {xi) = 9.0, and the angle change is A6 = 109.3°. 
Note also that the direction of J z is reversed. 





a) 

nT 




^ XSSS '\ ********** 
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b) 
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FIG. 2. Two examples of pair-plasma equilibrium state in a 
1-D configuration A+ = A_ = 1. There are components of 
B and J in both directions perpendicular to x, with a finite 
transition width dictated by the equilibrium equations and 
the adiabatic condition. In (a), the magnetic field vanishes at 
the origin, while it retains a finite value in (b). 



1. Sheet-Pinch 

With a sheet-pinch type equilibrium in mind, take V x 
B = aB, with a a scalar function of x. Then inserting 
P + = P = P into equation (|35|) shows that 

P = e B. 

2n 

For this case, B 2 and P 2 are both constant, which means 
that n/G and a are both constants as well. The 1-D so- 
lution is the sheet pinch B = B (sm(ax)y + cos(ax)z). 
However, in order to maintain P and B parallel, we re- 
quire that 



n ( aG 

2 g 



1 _ n f aG\ 







,G 



e 2 a 2 Xa 



1 = 0. 



Thus the temperature and density determine the scale of 
the field. If e is small (skin depth << system size), we can 
take a — 1/A, none of the relativistic effects matter, and 
the A scale is set by the system boundary conditions. The 
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momentum is small compared with the magnetic field 
strength in this limit. 

However, if skin depth effects are important, 

Real solutions imply that A > 2eyjG/n. The fact that 
one solution diverges as e — > shows the singular pertur- 
bation nature of the underlying equations. 

V. DISCUSSION 

We have shown that the evolution equations for 
charged, relativistic, homentropic fluids can be expressed 
in a vortex form, using a minimal coupling method hav- 
ing the same form as the canonical momentum of single 
particle motion, but with the inclusion of thermal effects. 
This places the magnetic field and fluid vorticity on the 
same footing. Following this formalism, we found that 
the effective Hamiltonian constructed from this canoni- 
cal momentum leads to a well-posed variational principle, 
which allowed us to find relaxed equilibrium states of the 
relativistic equations of motion. These equilibria can ex- 
hibit structures on multiple spatial scales, and could be 
very useful in the modeling of hot astrophysical flows. 
We have made no claims as to the stability of the states- 
indeed we expect the configurations containing counter- 
streaming species to be unstable in the kinetic regime. 
However, such states could be useful as starting points 
for studies of relativistic magnetic reconnection. 

The class of one dimensional solutions has a wide vari- 
ety of possible configurations depending on the symme- 
tries between the fluids. The multitude of these states 
must be narrowed down based on the boundary condi- 
tions particular to a given application. One possible ap- 
plication of this type of configuration is the modeling 
of the 'striped wind' that results in Pulsar Wind Neb- 
ulae (PWNe) with unaligned rotating dipoles^ ' 22 i 23 . In 
this configuration, an expanding entropy wave consists of 
low entropy regions of constant magnetic field separated 
by high entropy, neutral pair-plasma current sheets. A 
typical PWNe of this type can have 7G ~ 100, mak- 
ing a relativistic treatment necessary. It is expected that 
the dissipation within the magnetic field transition region 
happens on a much shorter timescale than the expansion 
of the wave, enabling the application of equilibrium the- 
ory. The reversed field configuration described in Section 
IIV Al provides a detailed and physics-rich description of 
this transition region, evaluated in the frame moving with 
the entropy wave. Although we have presented here solu- 
tions with constant entropy inside the hot, dense transi- 
tion region, this assumption can easily be relaxed in the 
1-D case. Outside of the transition regions, the constant 
field portion of the wave can be described by the case 
presented in subsection IIV A II and the entropy can be 



much lower than in the transition region. The magnetic 
pressure in the constant field regions balance the thermal 
pressure within the current sheets, as in previous studies. 
We will further examine the application of the formalism 
presented in this paper to these configurations in future 
work. 



Appendix A: Conserved Quantities 

Below, we show under what conditions the variational 
constraints used in the minimization principle are con- 
served. We use the general conservation principle that 
if a 4- vector V" is divergence free (d u V v = 0), then the 
spatial integral of the time component is conserved, pro- 
vided surface terms vanish. The total energy is automat- 
ically an ideally conserved quantity because the closure 
J" = E s ^ implies that d^Z = 0. 

To define the total helicity of each species, we take 
the dual of M, MP" = l/2e a0 > lu M a p, and construct the 
helicity 4-current 

The divergence of this 4-current is 

= -2fi- (E + c/gQ), (AI) 

where we have used the fact that the divergence is a 
scalar to evaluate in the lab frame, and we use the same 
definitions as in Section ITT] From (fTTj) . we have the lab- 
frame value: 

T 

(E + c/qQ) = u x n Vcr, 

97 

which allows us to write 

bjt* = 2— (O- VW. (A2) 
97 

Thus we see that helicity will be conserved if the entropy 
is constant on lines of total vorticity, which is true for 
isentropic plasmas in the case that u||fi. In the limit 
of small mass, this implies that entropy is a magnetic 
flux function, a standard result for MHD. In addition, 
we note that if there exists an equation of state such 
that the thermodynamic term can be expressed as a full 
gradient, 

(T/ 7 )W = VC, (A3) 

then due to the fact that V • ft = 0, the 4-divergence of 
the helicity can still be written as a total divergence 

d^K" = • vc = v • Qnc) • ( A4 ) 

Then the species helicity 
h= [ K°d 3 x= f[B+(c/q)R]-[A+(c/q)P}d 3 x, (A5) 
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is a conserved quantity. Note that the equation of state 
condition (|A3p is identical to the conservation of vortic- 
ity condition in equation (|12|) ; namely that the baroclinic 
term vanishes. Thus we conclude that whenever the flow 
preserves total species vorticity J7, the helicity h is con- 
served. Note that in the 1-D examples of section [IV] the 
baroclinic term always vanishes regardless of the form of 
entropy. 

This helicity contains a term involving the momentum 
times the curl of the momentum (~ v 2 k). Thus it has a 
higher order of spatial derivative than the total energy, 
which contains only the momentum squared (~ v 2 ). It 
is for this reason that the helicity is more fragile than 
energy to dissipation, and we must use the variational 
principle (fl4l) . 
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